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Abstract 

We present a generalization of the representation in plane waves of Dirac delta, 5(x) = 
~Tk I^oo e~ lkx dk, namely 6(x) = 4^ f^°oa e q usm § the nonextensive-statistical-mechanics q- 

exponential function, e tx = [1 + (1 — q)ix] 1 ^ 1 ~" ^ with ef = e ix , being x any real number, for real 
values of q within the interval [1,2[. Concomitantly, with the development of these new represen- 
tations of Dirac delta, we also present two new families of representations of the transcendental 
number tt. Incidentally, we remark that the g-plane wave form which emerges, namely, e lkx , is 
normalizable for 1 < q < 3, in contrast with the standard one, e lkx , which is not. 



1 



I. INTRODUCTION 



Dirac delta is a distribution that is used in almost all branches of physics. Various 
representations of it have been discovered along the time. For example, it can be represented 
as a limit of a Gaussian or as a linear combination of plane waves, being the last one strongly 
related to the Fourier transform (FT), as we will show later. 

Dirac delta, S(x), obeys the following fundamental property: 

f{x)S{x)dx = f(0), (1) 
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where / : K — > C is a well-behaved function. From the equation above, we can see that if 
f(x) = 1 Vx £ M, we get the normalization condition 



5{x) dx = l. (2) 

J — oo 

Also, choosing f(x) = f(0)e lkx in ([T]), we obtain 

5{x)e ikx dx = l, (3) 



i.e., the FT of 5(x) equals one. Therefore, using the expression of the inverse FT we obtain 
the following representation of Dirac delta: 

1 f°° 

S(x) = — / e~ ikx dk , (4) 
2vr J-oo 

which can be interpreted as a linear combination of plane waves. We can rewrite the ex- 
pression above as 



1 ' L 



S( x ) = — lim / e~ ikx dk , (5) 

Z7T L^oo J _ L 



then, Dirac delta also can be represented as the following improper limit: 

r/ x ,. sin(Lx) 

Six) = lim — i — '- . 6 

L^oo TXX 

In 1988, a possible generalization of Boltzmann-Gibbs statistical mechanics was proposed[l|. 
This new theory, sometimes referred to as nonextensive statistical mechanics has been 
satisfactorily applied to handle a large number of physical phenomena (usually, metastable 
or quasi-stationary states of systems that are not consistent with the ergodic hypothesis; 
for example, systems in which long-range interactions or strong-correlations exist) 



Furthermore, the elaboration of nonextensive statistical mechanics required the general- 
ization of some mathematical functions (exponential, logarithm, etc.), operators (sum, 



product, Fourier transform, etc.) and theorems (central limit theorem) 18]. Particularly, 
the generalization of the exponential function, namely, the g-exponential function is defined 
by 

ej = [1 + (1 - q)x)t« (ej = e*) , (7) 

for any where the symbol [y} + means that [y} + — y, if y > 0, and [y] + — if y < 0. 

For pure imaginary ix, e l q x can be defined to be the principal value of 

ef=[l + (l-q)ix]£i (ef = e-). (8) 

The main purpose of the present paper is to generalize the representation in plane waves 
of Dirac delta, introduced in equation (j2J), using the g-exponential function defined above. 



II. REPRESENTATION OF DIRAC DELTA IN g-PLANE WAVES 
A. Proposition 

Let us introduce the following quantity: 

1 f°° 

S q {x) = ^-y J e"^ d£ with q e [1, 2[ , (9) 

which can be interpreted as a linear combination of g-plane waves, where c(q) is a constant 
that may depend on q. We intend to show later that S q (x) = 5(x) for all 1 < q < 2. 
Analogously to (jHJ), we may write 

5 q (x) = -L lim / e"** with q e [1, 2[ , (10) 
c(g) A->ocJ_ A 

therefore, by integrating, we can represent 5 q (x) as the following improper limit: 

9 sin < 2=2 arctan[(g — l)Ax] [ 
W = 7^ rrrlim Lr i=zT~ withge]l,2[. (11) 
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B. The normalization constant l/c(g) and the transcendental number tt 



The constant c(q) must be equal to 2tc at the limit q — > 1 + . Furthermore, c(q) can be 
found from the normalization condition (J2]). Thus, we have 

2 poo sin < arctan[(g — I) Ax] > 

(2-9) A ^°°i-oo x[l + (g-l) 2 A 2 x 2 ]^y 
Using the change of variables z = (q — I) Ax we obtain 

o poo sin ( arctan z ) 
C ( q ) = ^—H m / \1 r-^dz. (13) 

(2 - q) A^oo J_ 00 ^ + ^2)^17 

As the integral does not depend on A, the limit symbol can be omitted. Therefore, we can 
write 

2 r°° sin [2a (q) arctan z\ 
(2~^q) J ^ Z (l + Z 2y 



c (<?) = / — . dz > ( 14 ) 



where 



a{q) = 2^T) • ( 15 ) 

We easily verify that a :]1, 2[c M — >■ M + is a monotonically decreasing function of q. 

In order to solve analytically the integral in f fl4|) . let us restrict to integer or half- integer 
values for a(q), more precisely, 1/2, 1, 3/2, . . .. This implies that q will be allowed to assume 
just certain rational values within the interval ]1, 2[, namely q = 3/2, 4/3, 5/4, . . .. Using the 
change of variables z = tan 8 in Eq. ()14p . we obtain 



= ^_ r<> s m [2 aM S](cosS)^>- 
By using now the relation ()A.4j) proved in the Appendix, the expression above yields 



c(g) 



2-g 

* fc=0 



V (-l) fc f 2a(g) ) f ^(cos^) 4 ^- 2fc - 2 (sin^) 2fc . (17) 



We remind that the Beta function, B(x,y), is defined by 

r-7r/2 



B{x,y) = d<f)2{cos<f)) 2x ~ 1 {sm<f)) 2y - 1 with x > and y> 0, (18) 
which is related with Gamma function by 

«*»-W$- (l9) 



Therefore, using the expressions of Beta function shown above, Eq. (|T7|) can be written as 

c[q > 2-q £j 1 ' r(2i + 2)r(2a(?) - 2k) ' { ' 

Let us rewrite now the expression above as 

<q) = ^—S nq , (21) 
z q 

where 

I 22±i|-l 

L 2 J V (n — k — -\V (k -i- -\ 

Sn =n q y (-l) fc V, * 2 / 1 T, 2j with n g = 2a(g) G N . (22) 
q ^ T(2k + 2)T(n q -2k) q w v ; 

When n q = 1 (which corresponds to g = 3/2) we obtain straightforwardly that 5*1 = ir. Also, 

it is straightforward to verify that S 2 , S 3 , £4 are equal to 7r. Using a symbolic computation 

software, we also verified that from n q = 1 to n q = 5000 (q = 5002/5001), S n = ir. Hence 

we state the following hypothesis: 

W" 1 r(n-fc-i)r(A:+i) 

7r = n V -1 fe — i— 2 ^ V 2y VnGN. 23 

^ v 7 T(2k + 2)r(n - 2fc) v ; 

We thus found a countable infinite family of representations of the transcendental number 



7r (see also 



19|). 



Using relation ( 123"]) in f l21|) . the expression of c(g) becomes 

= • (24) 

z q 

In addition to the above, this relation has been checked numerically to be correct not only 

for certain rational values of q within the interval [1, 2[, but for all real numbers within that 

interval (see Fig. [1]). Therefore, we conjecture that the integral which appears in (I14p equals 

7r for any value of q within that interval. Consistently, we obtain another infinite family of 

representations of the number ir, namely 

sin(2r arctan z) 
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dz Vr6l + . (25) 



This family is non countable and contains Eq. (T2"3"|) as a particular case. 

Finally, expressions (Q and ffTTj) of 6 q (x) become respectively 

2-q 



and 



/oo 
e-^di with q e [1,2 [ (26) 
-00 



sin < arctan[(g — l)Ax] > 

6Jx) = lim l - ^—L with q G]l,2[. (27) 

A ^°° nx[l + (q-l) 2 A 2 x 2 ]^Tv 
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FIG. 1. The blue dots were numerically obtained using expression (|14p . whereas the red continuous 
curve is the plot of c(q) given by equation §R 



C. Dirac delta behavior of the distribution 5„(x) 



Let us define the following distribution 



sin < ^1 arctan[(g — l)Ax] [ 
A q (x,A) = ^— 2 _ q J withge]l,2[, (28) 



irx [1 + (q-l) 2 A 2 x 2 ]^v 
which is related to S q (x) through 

SJx) = lim AJx,A) . 

A^oo 



(29) 



The plot of such a distribution (see Fig. |2]) indicates that in the limit A — > oo, A q (x, A) will 
present a divergence at the origin and will be zero for all i ^ 0, i.e., at first glance, 5 q (x) 
appears to be a representation of Dirac delta. 

Let us now consider an analytic function, / : dom/ C R — > C, which can be expanded 
in Taylor series around the origin such that the expression 



k=0 



m = E 

is valid for all x £ dom /. Then we have 

f(x)S q (x) dx 



/ (fc) (Q) fc 

k\ 



f(x)S q (x) dx. 



(30) 



(31) 



dom / 



Replacing f(x) by its Taylor series, this expression yields 



/W(0) 



x 



_1 sin ||zf arctan[(g — l)Ax]| 



/ f(x)5 q (x) dx = lim V , 

Jdomf A ^°°t^0 kl J <iomf TT [1 + (q - l) 2 A 2 X 2 }^^ 



dx, (32) 
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FIG. 2. Plot of A 3 / 2 (x, A) for different values of A. Similar results are obtained for any value of 
g€]l,2[. 



in which we must remark that q belongs to the interval ]1, 2[. If dom / is a bounded interval 
of R, i.e., dom / —]a, b[, with a < b, then using the change of variables z = (q — l)Ax, we 
obtain 

'(g-i)A6 z k ~ x sin ^^5f arctan , 
^^ (q-l) k A k k\J (q _ 1)Aa ~ 



[ f(x)8 q (x) dx = lim V - ^ lU °, fc ,, / 
Jdomf A ^°° t^Al- ^> k A k k\ J u 



The first term of the sum that appears above is 



7r (1 + 2 2 ) 2 ( 2 "-D 



dz. (33) 



r(q-i)Ab sin ( ^| arctan z 
/(O)lim / — ^ > 

A ^°°J( 9 -l)Aa 7TZ(1 + 2 2 ) 5 <^ T ) 



(34) 



If0<a<6ora<6<0, we straightforwardly see that expression above is equal to zero. If 
a < < b, then using relation (1251) we obtain that expression fl3ll) is equal to /(0). Finally, 
if we have either a = or b = (with a < b), then, also using relation (1251) we obtain that 
expression (13"1|) is equal to /(0)/2. 

In order to analyze the next terms of the sum given in (133]) . let us first rewrite them as 



lim / (fc) (0) J fe (A) with fc G N and fc > 1 , 

A— »oo 



(35) 



where 

Jk{q,A) 



(g-l)Afe Z 



fc-1 



sin f |zf arctan 



2-9 



dz with fc G N and k > 1 . (36) 



(q - l) k A k k\ J [q ^ 1)Aa ^(1+^2)2(^1) 

Jk{q,A) is a rapidly decreasing function of k (see Fig. [3]), which makes the sum given in 
( I3"3"j) converge, consistently with the finiteness of the domain of / in integral in Eq. f[3"Tj) . 



Moreover, from Fig. EJ we can infer that, in the limit A — > oo, Jk{q,A) — > 0. Therefore, Eq. 
implies 

' /(0) ifa<0<6 
/(0) /2 if either a = 0<6ora<0 = 6 (37) 
ifO0]a,6[. 



f(x)S q (x) d,X = < 
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FIG. 3. The /c-dependence of Jk(q,A), numerically obtained, considering b = —a = 1, for q = 1.4, 
and different values of A. From top to bottom, A = 10, A = 10 10 , A = 10 20 . 

In the case when dom/ is unbounded, i.e. if dom/ =]a, oo[, or dom/ =] — oo,b[, or 
dom/ = R, a similar analysis yields once again relation fl3T|) . Moreover, we numerically 
tested the validity of the mentioned relation using some types of functions and distributions 
(for example the Gaussian and the Lorentzian). Hence it seems reasonable to conjecture 
that, for a wide class of functions, S q (x) indeed is a representation of Dirac delta. Thus, we 
can finally write 

Z ~ 1 I „-ii{x-x') 



8(x — x) 



2tt 



e,-~ ~>de (?e[l,2[). 



(38) 



III. SQUARE INTEGRABILITY OF g-PLANE WAVES 



Let us consider the following function: 

ty(x) = Nef x = N [cos q {^x) + i sm q (£x)] with q e]l, 3[ and ^ , (39) 

which can be interpreted as a stationary g-plane wave, where the g-generalized trigonometric 
functions are defined, for any x G R, by (see also [20]): 

cos arctan[(g — l)x] | 



cosg x = Re (e q x ) 



[l + (q- l) 2 x 



21 2(g-l) 



(40) 



and 



sin, x 



Im ( e ») 



sin l^rj arctan[(g — l)x] 



[1 + (q - l) 2 x 2 } 



2(9-1) 



We illustrate these functions in Fig. HI 



cos 
1.0 




(41) 



sini ix 
l.Or' n 




FIG. 4. Top: cosi.i x (continuous curve) and cosx (dashed curve). Bottom: sini.i cc (continuous) 
and since (dashed). For 1 < q < 3, cos g z (sin g a;) is an even (odd) function of x. For 1 < q < 3, 
both functions cos q x and sin q x quickly decay when \x\ — > oo, in contrast with cosx and sinx. 



We will determine now the value of the constant N using the normalization condition 
given by 



ty*(x)^(x) dx = 1 . 



(42) 



Thus, we have 



-y roc 
AT2 = 



e q l &ef x dx 



(43) 



9 



which, using the definition of the g-exponential function given in (jSJ), can be written as 



1 r 00 1 

N2 Loo [1 + (a - 1) 



— dx . 



Using the change of variables tan# = (q — l)\£\x, this relation yields 

i i r /2 4-2, 

Therefore, we obtain that the normalization constant is given by 



N 



(44) 



(45) 



(?-i)l*|r A 



3-g 



(46) 



Let us emphasize that the function if)(x) = (plane wave) cannot be normalized, whereas 
g-plane waves, with q e]l, 3[, have a finite norm. 

0.6, 



N 



V~MT 




FIG. 5. The normalization constant N as a function of q. N goes to zero in the q — >■ 1 limit, thus 
recovering the well known non-normalizability of the plane waves. 



IV. CONCLUSIONS 

From the analytical and numerical results shown in section II, we conjecture Eq. (138]) . i.e., 
that S g (x) is indeed a generalization of the standard representation of Dirac delta in plane 
waves. Further research is welcome in order to establish which precise class of functions 
satisfy the relation (1571) . 

Concomitantly, we found two new families of representations, namely expressions (f2"5|) 
and (1251) . of the transcendental number it. We tested the validity of such expressions for 
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a set of values n G N and a G R. A demonstration is still required in order to formally 
establish these new families of representations of ir. 

A generalization of FT, namely, the so-called g-Fourier transform (g-FT) was developed in 
order to generalize the central limit theorem. The possible analytic expression of the inverse 
g-FT remains to be found. It is known that, using the representation in plane waves of Dirac 
delta together with the expression of the direct FT, it is possible to find the expression of the 
inverse FT. Consequently, we suppose that the present (/-generalization of the representation 
in plane waves of Dirac delta might be helpful in searching for an analytic expression of the 
inverse q-FT. Moreover, the present new representations of Dirac delta could be useful to 
handle some integrals that may appear in the analysis of certain physical phenomena. 

Finally, we prove a physically appealing property, namely that the g-plane wave form 
e l g X is square-integrable (in other words, normalizable) for 1 < q < 3, in contrast with the 
standard form, e lkx , which is not. 
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Appendix: Trigonometric identity 

We establish here an expression for sin[2o:(g)#], with 2a(q) G N and 9 G R, written in 
terms of sin6> and cos#. Firstly, we have 



sin[2a:(g)#] = Im (cos 6 + i sin 6) 



2a(q) 



(A.l) 



then, using binomial expansion we have 




(A.3) 



(A.2) 



(odd) 
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Finally, this expression can be rewritten as follows: 

sm[2a(q)9] = £ ° 19 (cos^^-^^sin^)* 1 (A.4) 

fc=0 V 7 

where we have used the floor function [_ J , defined, for any real number x, by [x\ = n such 
that n < x < n + 1, with n G Z. 
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